Investigation of the triple-alpha reaction in a full three-body approach 
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Background: The triple-alpha reaction is the key to our understanding about the nucleosynthesis 
and the observed abundance of 12 C in stars. The theory of this process is well established at 
high temperatures but rather ambiguous in the low temperature regime where measurements are 
impossible. Purpose: Develop a new three-body method, which tackles properly the scattering 
boundary condition for three charged particles and takes into account both the resonant and the 
non-resonant reaction mechanisms on the same footing, to compute the triple-alpha reaction rate 
at low temperatures. Methods: We combine the R-matrix expansion, the R-matrix propagation 
method, and the screening technique in the hyperspherical harmonics basis. Results: Both the 2f 
bound state and the 0^ resonant state in 12 C are well reproduced. We also study the cluster structure 
of these states. We calculate the triple-alpha reaction rate for T = 0.01 — 0.1 GK. Conclusions: 
We obtain the same rate as NACRE for temperatures above 0.07 GK, but the new rate is largely 
enhanced at lower temperatures (w 10 12 at 0.02 GK). The differences are caused by the direct 
capture contribution to the reaction when three alpha particles can not reach the resonant energies. 
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I. INTRODUCTION 

The formation of 12 C through the triple- alpha reaction 
was first suggested by Salpeter to explain the synthesis of 
heavy nuclei since stable isotopes of mass number 5 and 
8 can not be found in nature [1]. However, this reaction 
rate was too small to account for the observed abundance 
of 12 C in stars. Later on, Fred Hoylc predicted the ex- 
istence of a + resonance near the three-alpha threshold 
which significantly increased the reaction rate [2]. This 
resonant state and its properties were confirmed by ex- 
perimentalists [3, 4] shortly thereafter. Then, the two- 
step (sequential) mechanism of the triple-alpha capture, 

a + a -> 8 Be(0^) (1) 
8 Be(0+) + a^ 12 C(0+), (2) 

was accepted to explain the mystery of nucleosynthesis 
and the abundance of 12 C in helium burning stars. 

The first evaluation of the sequential (resonant) triple- 
alpha reaction rate was credited to Caughlan and Fowler 
[5] in which only the Hoyle resonance was considered and 
a Brcit-Wigner shape for the resonant cross sections of 
a — a and a— 8 Be was used. Ten years later, this reaction 
rate was improved by Angulo [6] when taking into ac- 
count the contribution of the 2^" resonance. This rate is 
included in NACRE: Nuclear Astrophysics Compilation 
of REaction Rates. At low temperatures, three alpha 
particles do not have access to the intermediate reso- 
nances, so one expects the three-alpha direct capture to 
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be dominant. The non- resonant triple-alpha rate tabu- 
lated in NACRE is obtained by a simple extrapolation 
of the sequential model to low energies. Since at low 
temperatures the measurements are impossible, a well 
founded theory for the three-alpha direct (non-resonant) 
capture is crucial. 

Ogata et al. [7] was the first group who attempted 
to tackle this problem without using an ambiguous ex- 
trapolation. In their studies, the continuum discretizcd 
coupled channel (CDCC) method [8] was employed to 
solve a three-body scattering equation, offering an op- 
portunity for the inclusion of both resonant and non- 
resonant mechanisms in the triple-alpha reaction. How- 
ever, the CDCC results [7] drastically differ from NACRE 
at low temperatures (20 orders of magnitude at T = 0.02 
GK) which significantly affects astrophysical studies and 
produces results that are incompatible with observation 
[9-13]. Applying the CDCC method to the triple-alpha 
problem can be challenging because the CDCC wavefunc- 
tion which is usually truncated can not reproduce the 
correct scattering behavior of three charged particles in 
the asymptotic region [14]. 

Another attempt using the hyperspherical adiabatic 
expansion method to solve the triple-alpha problem was 
made by the Madrid- Aarhus collaboration [15, 16]. How- 
ever, they were unable to perform the numerical calcu- 
lation at temperatures below T w 0.1 GK. Therefore, in 
order to estimate the non-resonant triple-alpha rate at 
low temperatures, they extrapolated a three-body Brcit- 
Wigner cross section for the three-alpha capture to low 
energies [17] in a similar manner as done in the sequential 
process [6]. The reaction rate from this work (named as 
BW(3B)) shows an increase of 7 orders of magnitude at 
T = 0.02 GK compared to NACRE. There are large dis- 



crepancies in the results of [6], [7] and [17], demonstrat- 
ing uncertainties and ambiguities in our understanding 
of the low temperature triple-alpha reaction. The aim of 
our work is to resolve this problem. 

At low temperature, the triple-alpha reaction proceeds 
through the + continuum states to the 2+ bound state 
of 12 C. In order to construct a reliable non- resonant reac- 
tion rate, we need a good description of both 12 C bound 
and continuum states. Recent microscopic theories such 
as the no-core shell model [18] and the Green's function 
Monte Carlo method [19] failed to reproduce the Hoyle 
state of 12 C by solving a true 12-body problem. Even 
though the fermionic molecular dynamics method pro- 
duces both the 2+ bound state as well as the 0^" reso- 
nance with significant triple-alpha configurations [20], a 
microscopic description for the non-resonant continuum 
states is currently not available. Given the dominant 
triple-alpha structure in both 12 C(2+) and 12 C(0^), it is 
reasonable to construct this as a three-body problem as 
done in [7, 15, 16]. 

Theories for the three charged particle problem have 
made significant progress over the last few decades. An 
accurate solution for a three-body bound state system 
can now be obtained (e.g [21]). However, there are still 
difficulties remaining for a three charged particle scatter- 
ing problem because no standard boundary condition for 
this system exists when considering the long-range effects 
of the Coulomb interaction. Some efforts have been made 
to solve this scattering problem but the results are lim- 
ited to cases where only one pair of particles have charge 
[22, 23]. Our triple-alpha problem is more difficult since 
all three particles are charged. In this work we employ 
the hyperspherical harmonics (HH) method in combina- 
tion with the R-matrix expansion and the R-matrix prop- 
agation technique to provide a good description of the 
a + a + a system at low relative energies with the cor- 
rect Coulomb asymptotic behavior, thereby enabling us 
to compute the reaction rate at very low temperatures 
without extrapolation. Our new three-body method is 
named here as HHR. This article is an extended version 
of the letter in [24], in which more details on the theo- 
retical methods are presented and a more comprehensive 
analysis of the results is performed. 

We have organized this article in the following way: 
detailed theory on the method is presented in Sec. II; the 
main results are shown in Sec. III. In Sec. IV we discuss 
in depth the reaction dynamics and the conclusions are 
drawn in Sec. V. 



II. THEORY 

A. Hyperspherical formalism 

The hyperspherical harmonics (HH) method originated 
in atomic and molecular physics [25, 26]. Later, it was ex- 
tended to few-body systems in nuclear physics by Delves 
[27, 28]. The theory of the HH method was well de- 




■ • 

1 1 



FIG. 1: (Color online) Coordinates for three alpha particles 
in Jacobi set i. 

veloped for Borromean systems [29-31]. Borromean nu- 
clei are defined as systems of three particles which are 
loosely bound and have no bound states in any of the 
two-body subsystems. 11 Li and 6 He which have two neu- 
trons weakly coupled to the core, are typical examples of 
Borromean nuclei. In this work, we formulate our prob- 
lem as a Borromean system of three alpha particles. 

Let us consider a system of three nuclei with masses 
rrii. For each Jacobi set i, we define the scaled coor- 
dinates x, = ^^0- k r, and y t = yj ^+a^+aI R * ( see 
Fig. 1). Here, = a^ — &j is the relative radius from par- 
ticle j to particle k and R; the radius from 
particle i to the center of mass of the two-body subsystem 
(jk). The ratios Aj = rrii/m are dimensionless and the 
scaled mass m is usually taken as the nucleon mass. The 
hyperspherical coordinates are then defined as functions 
of Xi and j/j 

2 2,2 

P = x i + Vi , 

9i = arctan — . (3) 

Vi 

The hyper-radius p is invariant under translations, rota- 
tions, and permutations while the hyper-angle 9i depends 
on the selected Jacobi set. For the triple-alpha problem 
in which three particles are identical, all the Jacobi co- 
ordinate sets are equivalent, so we choose to work with 
Jacobi set i — 3 and drop the index i in our equations 
from now on for convenience. 

Because our problem involves only spin-zero particles, 
all the equations in HH coordinates are written for this 
simpler form. The general case of non-zero spin particles 
is presented with details in [32, 33]. The hyperspherical 
expansion separates the radial and angular dependence 
of the three-body wavefunction 

*™ =p -i J2 XKU y (p)^K lv (9)[Yi^Y ly ] LM . (4) 

Kl^ly 

In Eq. (4) XKlj y {p) are the hyper-radial functions, solu- 
tions of the coupled channels equations (5) . The hyper- 
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angular functions ip^ y (6) are the eigensolutions of the 
hyper-angle equation and are defined by Jacobi polyno- 
mials [32]. This expansion introduces a new quantum 
number, the hyper-momentum K, an extended concept 
of angular momentum for a three-body system. l x and 
l y are the orbital angular momenta corresponding to the 
Jacobi coordinates x and y. 

Using the HH expansion above, we arrive at a set of 
coupled channels equations in the hyper-radius coordi- 
nate: 



2m 
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(5) 

where A = K + 3/2 and 7 = {K,l x ,l y }. The coupling 
potentials V~ ( y(p) are defined as the sum of three pair- 
wise interactions Vjk{p, 0) plus a three-body force V 3 b(p), 
integrated over all variables but p: 



V 7 y( P ) = (O 7 (0,x,y)| V jk + V 3b \ny(9,x,y)) , (6) 

k>] = l 

where Q 7 (0,x,y) = \Y ix <g> Y ly ] are the hyper- 

harmonic basis functions containing all the angular de- 
pendence in Eq. (4). In the HH representation, both 
the diagonal and the off-diagonal couplings decay slowly 
as Z^,/p for a system of three charged particles [14], 
demonstrating the difficulty of our probem in the asymp- 
totic region. Z°y is constant for any given channel 7 and 
7' and has a non-trivial expression: 



Zfy=Z 2 C 2 ' 
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where Z and A are the charge number and the atomic 
number of an alpha particle; Si^ is the Raynal-Revai 
coefficient [34]; i = 1, 2, 3 represent three Jacobi co- 
ordinates (xi,y ; ) and is the hyperharmonic basis 
function Q 7 (#j, Xj, yj) in Jacobi i. Each matrix element 
is integrated over all angular variables. 

For neutral Borromean systems, the three-body bound 
state wavefunction decays exponentially at large distance 
[35]. When introducing Coulomb interactions the prob- 
lem becomes more difficult because a simple analytic ex- 
pression for the asymptotic bound state wavefunction 
of three charged particles does not exist. However, for 
a well-bound system, imposing the boundary condition 
X-y(p — > 00) — > is sufficient for the numerical calcula- 
tion [32]. 

Solving Eq. (5) for positive energies E, to obtain con- 
tinuum states for a system of three charged particles, is 



numerically challenging since it requires an exact bound- 
ary condition. Fortunately, there are techniques (which 
will be discussed in detail later) that enable us to neglect 
the off-diagonal potential couplings at large distances. 
This results in analytic solutions for the asymptotic ver- 
sion of Eq. (5) corresponding to the regular F^(r)j, np) 
and irregular Ga(ij 7 , up) three-body Coulomb functions. 
We can thus write down the asymptotic three-body con- 
tinuum wavefunction for a charged system in hyper- 
radius: 



X ™« (p) 



H A (jy 7 , np)5 1 y - H^(r] 7 , np)S, 



77 



(8) 



= G ±iF are the Coulomb functions describing the 
outgoing and incoming spherical waves and S^y is the 
scattering matrix [33] . Due to the nature of the scattering 
boundary condition, each partial wave x 77 ' requires two 
subscripts 7 and 7' to describe the outgoing and incoming 
channels respectively. The hyper-radial momentum k is 
given as a function of the three-body kinetic energy k = 
y/2mE/h 2 , where H is the Planck constant. The variable 
?7 7 is an equivalent of the Sommcrfcld parameter in two- 
body systems but is more complex and hyper-momentum 
K dependent. ?y 7 is obtained from the parameter Z^ 
through a relationship i] 7 = mZ^/fPn. 



B. R-matrix expansion 

Wigner and Eiscnbud [36] were the first to introduce 
the R-matrix method into nuclear physics to study reso- 
nances in nuclear scattering. Later on, this method was 
developed and extended in detail by Lane and Thomas 
[37]. The general idea of this method is to use an or- 
thonormal basis expansion inside an R-matrix box. An 
R matrix is then constructed to match to the asymp- 
totic wavefunction outside the box [33]. Even though 
the R-matrix method was originally developed for a two- 
body scattering problem, it has been generalized for a 
three-body system in HH coordinates by I.J. Thompson 
el al. [38]. Typically the R-matrix expansion provides a 
method that is numerically more stable than direct inte- 
gration methods for solving coupled channels equations 
with strong repulsive couplings [33]. 

By imposing a constant logarithmic derivative /3 = 
dln<j>P(p)/dp at a fixed radius p m , the so called R-matrix 
radius, the eigensolutions </> 7 (p) of Eq. (5) will form an 
orthonormal basis over a finite range [0,p m ]. The scat- 
tering wavefunctions at energy E are expanded in terms 
of (jjP(p) at eigenenergies e p : 



X 77 ' 



2J ^77' 



(9) 



In our calculation, the logarithmic derivative j3 is cho- 
sen to be zero. This highly reduces the complexity of 
the R-matrix propagation method (Sec. II C), therefore 
increases the efficiency and accuracy of the numerical 
performance when propagation is needed to a very large 
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radius. All the equations presented in this paper con- 
sider j3 — while more general formulae can be found in 
[33] . The R matrix at energy E relates the wavefunction 
X 77 ' (p) and its derivative: 



X77' (p)=J2p R ~n" ( E ) x'yy (P) • 



(10) 



7" 



The constant logarithmic derivative enables us to 
calculate the R matrix i? 77 »(i?) from the known eigen- 
functions 4> p at radius p m : 
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(11) 



where P is the number of poles used in our calculation. If 
we substitute Eq. (9) and Eq. (11) into Eq. (10), the only 
missing quantities to calculate the scattering wavefunc- 
tion are the expansion coefficients A? , . As discussed in 
Sec. II A, Eq. (8) is the asymptotic form of a continuum 
wavefunction for three charged particles when the off- 
diagonal couplings can be neglected. For sufficiently large 
p m the wavefunction will be in the asymptotic regime. 
Inserting Eq. (9) and Eq. (11) in Eq. (10) and combining 
with Eq. (8) we obtain the unknown coefficients A? , : 

p 7 " 

(12) 

All functions in Eq. (12) are evaluated at the R-matrix 
boundary p m . The scattering matrix Sy/y and the R 
matrix are directly related [33]. 



C. R-matrix propagation 

The triple-alpha system is driven by strong, long-range 
pairwise Coulomb interactions. Therefore, the R-matrix 
radius p m has to be very large to fully capture the physics 
of the problem. Solving the coupled channels equations in 
a large R-matrix box will lead to numerical instability. In 
order to overcome this difficulty, we employ the R-matrix 
propagation technique originally developed by Light and 
Walker for the atom-molecule scattering problem [39, 40]. 
First, the set of coupled channels equations (5) is solved 
in a small R-matrix box with a radius of p m . Then the R 
matrix at p m is propagated to p a 3> p m where matching 
to known Coulomb functions can be performed. This 
propagation method is well known for its fast convergence 
and high numerical stability. 

We rewrite Eq. (5) as 



d 2 X 7 7'(p) 

dp 2 



= E*r 

7" 



X 7 "7'(P) 



(13) 



The interval from the R-matrix radius p m to the match- 
ing radius of the three-body asymptotic wavefunction p a 



is divided into sectors. We choose the size h p of sector 
p to be sufficiently small that the interaction within is 
considered to be constant A 7 (p) in each channel 7. We 
diagonalize V 7 y by solving the equation: 



(T* ) T V( Pp ) T p = A(p) 2 



(14) 



where p p is taken at the center of sector p. The wave- 
function and its derivative at the left and right boundary 
of sector p can be related through a local diagonal repre- 
sentation of the propagating functions G p : 



(15) 



The subscripts R, and L imply evaluations at the right 
and left side of the sector boundary, respectively. The 
propagating functions are expressed as Gi P = T p gi p T p , 
where g p is a simple function of A 7 (p): 
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The R matrix is then propagated from sector p— 1 to the 
next by: 



Pi 



j (G 2 p [G 4 p 



PR 



RP -i]-i Ga P _ Gi j 



Eq. (17) enables us to propagate the R matrix from a 
small radius p m to a distance p a large enough such that 
the asymptotic matching can be performed and hence 
the wavefunction at the R-matrix boundary can be cal- 
culated. Then the propagating functions G p are used 
again to reconstruct the continuum wavefunctions to the 
desired radius 



- 1 xT 1 



Xl = G-l p ~ [( 



.(18) 



In our calculations, the logarithmic derivative j3 is al- 
ways kept at zero so that Eq. (17) and Eq. (18) remains 
valid. 



D. Screening 

The presence of very narrow two-body and three-body 
resonances in addition to the strong, long-range Coulomb 
interaction in our problem creates numerical instabilities 
in the propagated three-body scattering wavefunctions. 
We overcome these difficulties by introducing a Woods- 
Saxon screening factor [1 + exp(( ( o — /5 S crcon)/ a scrocn)] 1 
in the off-diagonal potentials of Eq. (6). The screening 
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FIG. 2: (Color online) Our three-body method is divided into 
four steps in which we employ the R-matrix expansion, R- 
matrix propagation and screening technique in the Hyper- 
spherical Harmonics coordinates. 



radius p scrccn is chosen sufficiently large to gain conver- 
gence, yet small enough to ensure numerical stability. 

Fig. 2 summarizes the different steps in our method. 
We divide the hyper-radial space into four regions. The 
Hyperspherical Harmonics R-matrix expansion (HHR) 
is first applied in a small box of radius p m (~ 50 fm). 
We then use the R-matrix propagation method to much 
larger distances. All Coulomb couplings are included 
from p m to Pscreen (~ 800 fm) . After that, the off- 
diagonal couplings are screened up to p a (~ 3000 fm) 
where it is safe to perform the asymptotic matching. 



E. Three-body reaction rate 

For a given three-body radiative capture reaction a + 
b + c — > D + 7, the reaction rate R a bc at a given energy E 
is calculated through the photo-dissociation cross section 



h 3 

Rabc(E) = pi N\ G abc , D ^ 



8tt 



E 2 
7 r 7 (25 7 ) 



(19) 



The photon energy E y and the three-body kinetic energy 
E are related by E 1 = E + \Ejj\. Here Ejj indicates the 
bound state energy of the final nucleus D relative to the 
three incident-particle a, 6, c threshold. p a b and p a b,c 
are the reduced masses of the a + b and (ab) + c systems, 
respectively, p is the number of identical particles and 
Na is Avogadro's number. The statistical factor G a bc,D 
depends on the spins of the nuclei a, 6, c and D through 
the relationship G abc , D = {2 j a +iWJ?+iH2j c +i) ■ 

The energy averaged reaction rate (R a bc){T) (a func- 
tion of temperature) is an important quantity in astro- 
physics. It is obtained by integrating R a f, c (E) over the 
three-body Maxwell-Boltzmann distribution which dif- 
fers from the two-body case: 



(Rabc) = 



1 1 



2 {k B Tf 



RabcE 2 e ■ 



^ dE , (20) 



where kg is the Boltzmann constant. Throughout this 
paper, we will refer to the energy averaged reaction rate 



(R a bc)(T) as the triple-alpha rate. 

The photo-dissociation cross section a 1 in Eq. (19) 
depends on the electromagnetic transition strength 
dB(EX)/dE through which the reaction occurs: 



(2tt) 3 (A+1) (E. 
A[(2A+1)!!] 2 
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(21) 



At low energies, the triple-alpha reaction proceeds 
through a quadrupole transition from the + continuum 
to the 2+ bound state in 12 C. The formula in the hy- 
perspherical basis for the transition strength function 
d_B(£^2)/d£' is given in the Appendix. A detailed deriva- 
tion of the three-body quadrupole transition and reaction 
rate is presented in [41]. 



III. RESULTS 

A new code named HHR3a is generated to solve Eq. (5) 
for bound states and continuum states of the triple-alpha 
system. It is developed from the program FaCE [21] and 
STURMXX [42] which are originally designed by I.J. 
Thompson et al. for a core + n + n system. Our cur- 
rent problem, involving three identical alpha particles, 
has different symmetries and is more difficult because of 
the long-range Coulomb interaction. An alpha particle is 
considered as a boson of spin zero, thus the wavefunction 
must be unchanged for any permutation of these two par- 
ticles. Only partial waves with l x = even will contribute 
to the total wavefunction. Here, l x is the relative angu- 
lar momentum of the two identical alpha particles being 
interchanged. The new code HHR3a is corrected for the 
symmetry properties introduced by the three identical 
bosons. In addition, the correct Coulomb asymptotic 
wavefunctions are used instead of plane waves in STUR- 
MXX [42] since our problem involves charged particles. 
We implement the screening technique in this new ver- 
sion of the code and considerable modifications are also 
introduced in the R-matrix propagation method. 



A. Interactions 

In our study, the triple- alpha problem is reduced 
to an effective three-body problem in which the two- 
body alpha-alpha interactions are phenomenological and 
a three-body force is needed to account for the fact that 
the alpha particle is not a fundamental particle. We use 
the Ali-Bodmcr potential for the alpha-alpha interaction 
[43] as modified by Fedorov et al. [44] : 



V aa = (125PU + 20P l=2 )e- 



,-r'/1.53 2 _ 30 .l 8e -r a /2.85'_ 

(22) 

The Z-dependent potential has a very strong repulsive 
core in the s-wave to simulate the Pauli exclusion prin- 
ciple in the a — a system. This interaction reproduces 
successfully the low energy phase shifts as well as the Of 
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resonant state of 8 Be at 0.093 MeV. 

A three-body force is defined by V^p) — Vae~ p / pa 
and added to the Hamiltonian to reproduce the exper- 
imental bound state energy of the 2\ state as well as 
the Hoyle resonant state 0^ in 12 C. We take po = 6 fm 
as it corresponds to the position where the three alpha 
particles touch each other. We adjust Vo to the value 
of -15.94 MeV (-19.46 MeV) to generate the 2+ bound 
state (the 0^" Hoyle resonance) at the correct binding en- 
ergy. Vo = —19.46 MeV is then used to calculate all the 
+ continuum states. 

In addition to the nuclear interaction we include the 
Coulomb potential: 



2„2 



V^ ul (r) = Z 2 e 



(s 2rL.) 



- r < r C oui 

r > r C oui ■ 
(23) 

Here, Z is the charge of an alpha particle, r is the distance 
between two alpha particles, and rcoui is the Coulomb 
radius which is taken as twice the alpha particle radius 
rcoui = 2.94 fm. 

In Table I, we summarize the parameters for the three- 
body interactions used in our calculations to reproduce 
the experimental binding energy for the 2+ bound state 
and the resonant energy for the Hoyle state 0^ of 12 C 
[45]. 



B. The 2J bound state 

A three-body bound state is obtained by solving the set 
of coupled channels equations (5) with a boundary condi- 
tion that the wavefunction goes to zero at large p. These 
coupled equations are solved by expanding the hyper- 
radial wavefunction in terms of the Laguerre basis [46] 
which forms a complete and orthogonal set. The matrix 
diagonalization is then applied to obtain eigenenergies 
and eigenvectors. We use the modified Ali-Bodmer in- 
teraction [43, 44] in combination with a three-body force 
as discussed in Sec. Ill A to obtain the 2+ bound state 
for 12 C in agreement with experiment. 

Fig. 3 and Fig. 4 show the convergence of the three- 
body binding energy and the rms radius of the 2|~ bound 
state of 12 C with the size of the model space represented 
here by the maximum hyper-momcntum K max . A given 
_ftT max determines how many channels are included in the 
wavefunction expansion Eq. (4), e.g., i^ max = 20 pro- 
duces 36 channels in the expansion. From Fig. 3, we 
see that the bound state energy starts converging at 



TABLE I: The three-body interactions to reproduce the ex- 
perimental energies [45] of the 2^ bound state and 0^" reso- 
nant state in 12 C 





Vb(MeV) 


Po(fm) 


E(MeV) 




-15.94 


6 


-2.875 


oi 


-19.46 


6 


0.380 




FIG. 3: Dependence of the three-body binding energy of 
the 12 C (2^~) bound state on the maximum hyper-momentum 
Jfmax which determines the size of the model space. Any 
K < 7f max will be included in the wavefunction expansion 
Eq. (4). 



2.47 



2.465 - 



2.46 



2.455 



2.45 




FIG. 4: Dependence of the r rms radius of the 12 C(2+) bound 
state on the maximum hyper-momentum if m ax. The physical 
meaning of A" m ax is described in Fig. 3. 



-Kmax = 12. Therefore, energy convergence is guaran- 
teed by choosing JT max = 20. The bound state energy 
of the 2f state of 12 C converges to the experimental 
value E = —2.875 MeV [45] with respect to the three- 
alpha threshold. Our final value for the rms radius is 
2.459 fm which is close to the result obtained in a micro- 
scopic calculation done by Chernykh et at (2.50 fm) [20]. 
This agreement validates the three-body approximation 
for the 2\ bound state. In support of this is also the fact 
that in [20], the triple-alpha configuration is dominant 
for this 12 C state. 
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FIG. 5: (Color online) The density distribution for the 2^ 
bound state of 12 C. The brighter color corresponds to the 
higher density distribution. 

We next construct the density distribution function 
from the 2 + bound state wavefunction to study its struc- 
ture: 

P(r,R) = J |*(r,R)| 2 drdR. (24) 

Here r is the radius between two alpha particles and R 
is the distance from their center of mass to the third al- 
pha. In Fig. 5 we present the density distribution P(r, R) 
for the 2^~ bound state of 12 C. The dominant configura- 
tion for this state is an equilateral triangle in which each 
pair of particles is ~ 3 fm apart. This configuration is 
represented in Fig. 9c. 

C. The 0+ Hoyle state 

The three alpha particle system is driven by a strong 
Coulomb interaction of which the off-diagonal couplings 
are long-range in the HH representation. Therefore, solv- 
ing the coupled channels equations (5) for positive ener- 
gies E is not a trivial task. The HHR method which com- 
bines the R-matrix expansion and the R-matrix propa- 
gation in the HH basis is employed to overcome this diffi- 
culty. We again use the modified alpha-alpha Ali-Bodmer 
potential [43, 44] and adjust the three-body force to re- 
produce the experimental Hoyle resonant state (details 
of interactions are shown in Sec. Ill A). 

Fig. 6 shows the convergence of the resonant energy 
of the 0j~ Hoyle state of 12 C with the maximum hyper- 
momentum -ftT max . As expected, we find that the conver- 
gence for the Hoyle state is slower than for the bound 
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FIG. 6: Dependence of the three-body energy of the 12 C(0j) 
Hoyle resonant state on the maximum hyper-momentum 




E(MeV) 

FIG. 7: The quadrupole transition strength d_B(_E2)/d_E as a 
function of the three-alpha kinetic energy E. 



state. The result begins to converge at K max — 26 and 
this is the point we choose to fit the three-body force to 
reproduce the experimental energy of 0.38 MeV [45] with 
respect to the three-alpha threshold. The uncertainty of 
the resonant energy is 4% as a result of comparing the 
value at if max = 26 and its extrapolation to K max — > oo. 

The + continuum states of 12 C are calculated to 
very low energies E. The quadrupole transition strength 
Eq. (Al) is then constructed using the 2+ bound state 
wavefunction and the + resonant and non-resonant con- 
tinuum states. Therefore, in our calculations we treat 
the resonant and non-resonant process on the same foot- 
ing. Fig. 7 plots the quadrupole transition strength as 
a function of the kinetic energy for the three interact- 
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FIG. 8: (Color online) The density distribution for the 0^" 
Hoyle resonant state of 12 C. The brighter color corresponds 
to the higher density distribution. 



ing alpha particles. The curve peaks at the measured 
Hoyle resonant energy {E = 0.38 MeV). As we expect, 
the strength function decreases as it approaches the lower 
energy regime. Around 0.05 MeV there is a sharp reduc- 
tion of the transition strength below which the formation 
of 12 C becomes unlikely. 

By using the wavefunction at the resonant energy we 
are able to construct a density distribution function as 
in Eq. (24). This quantity illustrates the spatial struc- 
ture of the Hoyle state, a topic of trememdous interest 
for many years [20, 47-50]. 

Fig. 8 depicts the density distribution in a color con- 
tour plot where the brighter color corresponds to a denser 
distribution. The dominant configuration for the Hoyle 
state is the prolate triangle in which two alpha particles 
are near each other (~ 3.7 fm) and further away from 
the third alpha (~ 5.2 fm), in agreement with the find- 
ings in [44, 50] (presented in Fig 9. a). We also observe 
two smaller maxima for the density distribution function. 
One of them supports an oblate triangle configuration of 
three alpha particles (Fig 9b) , where two alpha particles 
are ~ 7 fm apart and their center of mass is ~ 1.5 fm 
from the third alpha. The other small maximum indi- 
cates an almost equilateral triangle with a distance of 
~ 3 fm between any pair of alpha particles (Fig 9c). 
The weights of these two configurations are about half 
of that of the prolate triangle configuration. In conclu- 
sion, we find three different configurations in the Hoyle 
state: the prolate triangle is dominant, but the oblate 
triangle (almost chain-like) and the equilateral triangle 
also contribute significantly Fig 9 portrays accurately 
the configuration of three alpha clusters only in relative 
magnitudes of r and i?, and the average angle between 




(b) 

FIG. 9: (Color online) Different configuration of the triple- 
alpha system: the prolate triangle (a), the oblate triangle 
(b), and the equilateral triangle (c). 



these two vectors requires further investigation. 



D. Convergence and uncertainty 

It is important for us to carefully study the convergence 
properties of our method, given the well known difficulty 
of the problem. With the help of the R-matrix propaga- 
tion method we are able to perform calculations out to a 
very large radius and obtain stable results with Coulomb 
interactions only. When nuclear interactions are intro- 
duced, the very narrow resonances occur in both the two- 
body and three-body systems, and lead to a degradation 
of our propagation technique and increasing numerical 
instabilities, especially in the turning point region. As 
mentioned in the theory section, we tackled this prob- 
lem by screening the off-diagonal couplings by a Woods- 
Saxon multiplying factor [l+exp((p-p scrccn )/a scrccn )] _1 . 
However, it is necessary for us to ensure that important 
physics is not left out. 

In Fig. 10, we present three different calculations of 
the reaction rate as a function of temperature for a sys- 
tem of three alpha particles when only Coulomb interac- 
tions are included. For this case, there are no resonances 
in either the two-body or three-body system, thus the 
propagation technique is stable and we are able to de- 
termine the impact of introducing screening potentials. 
We perform a calculation using the screening technique 
with /O scrccn = 800 fm. This calculation fully converges 
at Kmax — 26 and is represented by the dashed line in 
Fig. 10. This result is then compared to two calculations 
with no screening at all. One is done with K max = 20 
and the other with i^ max = 26 to ensure generality and 
the convergence of the method. Uncertainty of the reac- 
tion rate due to the screening technique ranges from 35% 
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FIG. 10: (Color online) The reaction rate as a function of 
temperature for the fictitious three-alpha system in which 
only Coulomb interactions are introduced. Comparison be- 
tween screening and no screening calculations: with screen- 
ing and A* max = 26 (dashed), no screening and if max = 20 
(dot-dashed), no screening and if max = 26 (dotted). 



FIG. 12: (Color online) Convergence of the triple-alpha re- 
action rate at low temperatures (0.01-0.06 GK) with hyper- 
momentum -Kmax- Four calculations are presented: if max = 
20 (solid), if max = 22 (dashed), if ma x = 24 (dot-dashed), 
if max = 26 (dotted). The inset plot presents the triple-alpha 
rate for a large range of temperatures (0.01-1 GK) in which 
the contribution of the resonant triple-alpha capture occurs 
at T > 0.06 GK. 




FIG. 11: (Color online) Evaluations of the triple-alpha reac- 
tion rate using three different screening radius: p scrcC n = 600 
fin (solid), p S creen = 800 fm (dashed) , p S creen = 1000 fm (dot- 
dashed) . 
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FIG. 13: (Color online) Convergence of the quadrupole tran- 
sition strength &B(E2) / &E with hyper-momentum if max at 
E = 0.01 MeV. 



to 39% for temperatures within 0.02 — 1 GK. These are 
very small numbers given the many orders of magnitude 
involved in the problem. 

For our realistic problem of the triple-alpha where 
both nuclear and Coulomb interactions are included, we 
need to use the screening technique to obtain stable re- 
sults. The convergence study of our triple-alpha rate with 
varying p scro0 n an d if max is performed carefully. We use 
an R-matrix box size of 50 fm, with P = 50 poles in- 



cluded in the R-matrix expansion. Several values of the 
screening diffuseness are used and no dependence has 
been found, we thus fix the screening diffuseness to 10 
fm throughout this paper. The R matrix is then prop- 
agated out to 3000 fm, sufficiently large to perform the 
Coulomb asymptotic matching. Fig. 11 presents the re- 
sults using three different screening radii: 600, 800 and 
1000 fm. The triple-alpha reaction rate converges for 
Pscreen > 800 fm with an uncertainty of 5%. 



Fig. 12 shows the convergence of the triple-alpha re- 
action rate with hyper- momentum AT ma x- We perform 
four calculations of the rate as a function of temperature 
by varying K max from 20 to 26. -ft" max = 20 is not large 
enough to reproduce the experimental 0j Hoyle resonant 
state, therefore the reaction rate for this case is larger 
than other curves for T > 0.1 GK (see the inset plot). 
The triple-alpha rate is well converged at K max = 26. 

In Fig. 13 we look at the quadrupole strength func- 
tion, which is a principal ingredient for calculating the 
reaction rate, at a specific energy of E = 0.01 MeV. If 
dB(E2)/dE converges, then the cross section and the 
reaction rate will also converge. Since the dB(E2)/dE 
has an exponential behavior, namely y = A — Be~ Cx 
as a function of hyper-momentum K max , we are able to 
extrapolate to the converged value of dB(E2)/dE when 
K max —> oo. Even for a very low energy (E = 0.01 MeV), 
the extrapolated value of the rate differs from that at 
K max = 26 by only 4%. 

In order to study the sensitivity of the triple-alpha rate 
to the nuclear interaction and the three-body force, we 
consider employing the same Hamiltonian as [7] to calcu- 
late the + continuum. In [7] the alpha-alpha interaction 
is given as: 

V aa (r) = lOO.Oe- 1 - 2 / 1 002 - SO.Sbe^ 2 / 2 132 . (25) 

Unlike in [7], we find that a three-body interaction is 
needed to reproduce the relevant Hoyle state. As seen 
from Fig. 14, the new interaction produces a triple-alpha 
reaction rate at low temperatures 4 orders of magnitude 
higher than that obtained with the Ali-Bodmer interac- 
tion [43, 44]. While the Ali-Bodmer potential reproduces 
the a-a phase shifts, the interaction in [7] does not. It 
therefore provides an upper limit for the error associated 
with V aa ambiguities. 

Our final rate is normalized to the standard rate from 
NACRE for temperatures T > 0.5 GK, where the Hoyle 
resonance completely dominates. At this high tempera- 
ture, the triple-alpha reaction proceeds through the very 
narrow Hoyle resonance and its rate is proportional to a 
gamma decay width T. In NACRE, r = 3.7 x 10" 9 MeV 
is taken from experimental data [45]. We fit our calcu- 
lated cross sections to a Breit-Wigner shape and obtain 
a value of 2.86 x 10~ 9 MeV for the partial decay width. 
Therefore, the HHR rate is expected to be at least 1.3 
times smaller than the NACRE rate due to the different 
width. In practice, a factor of 2 is needed for normalizing 
to the NACRE rate. We consider this factor of 2 to be 
our conservative error counting for both the uncertainty 
of the decay width and the convergence of the problem. 

E. Rate 

With the 2+ bound state and the + continuum states 
of 12 C obtained above, we are able to compute the 
triple-alpha reaction rate at different temperatures. Our 
method enables us to treat the resonant and non-resonant 
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FIG. 14: (Color online) The sensitivity of the triple-alpha rate 
to the interactions: comparison between a— a Ali-Bodmer+3- 
body force (solid) and ct — ot interaction as in [7]+3-body force 
(dashed). 



processes on the same footing. We present our new triple- 
alpha rate in Table II for temperatures ranging from 0.01 
GK to 1.0 GK, beyond which we would need to consider 
other high lying resonances. 

The triple-alpha reaction proceeds primarily through 
either the resonant or the non-resonant path depending 
on the temperature of the stellar environment. In or- 
der to have a full description of the triple-alpha reaction 
mechanism, we estimate the energy range relevant for the 
reaction rate at a given temperature. These ranges of en- 
ergies contribute significantly to the integral of Eq. (20) . 
As shown in Fig. 15, the resonant energy at 0.38 MeV 
completely dominates the integration for T > 0.07 GK. 
For T = 0.06 — 0.07 GK, there is competition between 
the resonant and non-resonant processes. This marks the 
transition region between the two processes. The non- 



TABLE II: Our triple-alpha reaction rate (in cm 6 s 1 mol 2 ) 
after being normalized to NACRE for T > 0.5 GK. 
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FIG. 15: Range of energies relevant for the triple-alpha reac- 
tion rate at a given temperature. 



resonant capture mechanism dominates for T < 0.06 GK. 



FIG. 16: (Color online) Different evaluations of the triple- 
alpha reaction rate: comparing the Hyperspherical Harmonic 
R-matrix method (solid) with NACRE (dotted), CDCC 
(dashed) and the three-body Breit Wigner (dot-dashed). 



IV. DISCUSSION 

A. Comparison with other methods 

Fig. 16 presents our results (solid line) in comparison 
with other studies: NACRE, CDCC and BW(3B). The 
new rate resulting from the HHR method (after being 
rescaled as discussed in Sec. HID) agrees with NACRE 
(dotted line) for temperatures above 0.07 GK. Although 
the reaction rate is slightly reduced below 0.07 GK, it is 
significantly enhanced for T < 0.06 GK. We also obtained 
a very different temperature dependence in this region 
[24]. The results assuming an extrapolation of a three- 
body Breit- Wigner cross section to low energies BW(3B) 
(dot-dashed line) [17] have a similar behavior but the 
reaction rate increases to a lesser extent. The CDCC re- 
sults of [7] (long-dashed) largely enhance the triple-alpha 
rate for temperatures as large as 0.1 GK when compar- 
ing with NACRE. This effect is much stronger than what 
is seen in our studies. The fact that the HHR result 
agrees with NACRE for T > 0.07 GK leads to a negli- 
gible change in the evolution of stars around one solar 
mass [24]. This cannot be reproduced when using the 
CDCC rates. 

There exists a kink in the HHR curve around T w 0.06 
GK as seen in Fig. 16. This marks the transition be- 
tween the resonant and the non-resonant processes for 
which the temperature dependences are different. Above 
T w 0.06 GK the sequential (resonant) process domi- 
nates, while below there is mostly non-sequential (direct) 
capture. This agrees with the finding in Fig. 15. The cal- 
culations in [17] exhibits this same feature. 



B. Long-range Coulomb effects 

Since we obtain a large enhancement for the triple- 
alpha reaction rate at low temperatures, it is important 
to isolate the source of this effect. We therefore perform 
calculations of the 12 C(0 + ) continuum states for cases 
in which only the diagonal Coulomb couplings VfL (dot- 
ted) and the full Coulomb couplings VfL, (dot-dashed) 
are included. The results are then compared to calcu- 
lations including both nuclear and Coulomb interactions 
(V c + V N ) 17 > with the off-diagonal Coulomb couplings 
up to 30 fm (dashed) and 800 fm (solid). The reaction 
rate for each case is then constructed by using those scat- 
tering wavefunctions and fixing the 12 C(2+) bound state. 
Fig. 17 shows the results of these calculations. When only 
the diagonal Coulomb couplings are present, we are able 
to obtain an analytic solution of Eq. (5) (Sec. II A) . Our 
numerical calculation for this case agrees well with the 
analytic solution allowing us to test our implementation. 
The inclusion of the off-diagonal Coulomb couplings (dot- 
dashed curve) significantly increases the reaction rate at 
low temperatures. When both nuclear and Coulomb cou- 
plings are fully included in our calculation, we observe 
an increase in the reaction rate at high temperatures due 
to the resonant contribution. Comparision between the 
solid and the dashed curves confirms that off-diagonal 
Coulomb couplings drive the increase in the triple-alpha 
reaction rate at low temperatures. Fig. 17 indicates that 
the effect of off-diagonal long-range couplings are rela- 
tively small at high temperatures but very important in 
the low temperature regime. About 10 orders of magni- 
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FIG. 17: (Color online) The long-range Coulomb effects are 
shown in four different calculations: only diagonal Coulomb 
couplings (dotted), only Coulomb couplings (dot-dashed), 
both nuclear and Coulomb interactions with off-diagonal 
Coulomb couplings up to 30 fm (dashed) and a fully con- 
verged calculation with off-diagonal Coulomb couplings up to 
800 fm (solid). 



tude enhancement in the rate is found at T — 0.01 GK 
due to these effects, demonstrating the importance of in- 
cluding Coulomb correctly. 



C. Reaction dynamics 

In order to understand the mechanism for 12 C produc- 
tion in more detail, we rewrite Eq. (Al) as 



&B(E2) 
dE 



J fy<(r,R)drdR 



(26) 



where r is the radius between two alpha particles and 
R is the distance from their center of mass to the third 
alpha particle. 7* = {if 1 , Z* , Z* } indicates an incoming 
channel in the hyperspherical wave expansion for a scat- 
tering state. The spatial distribution of function / 7 < (r, R) 
at different three-body kinetic energies E contains infor- 
mation about the dynamics of the triple-alpha reaction. 
Fig. 18, 19, and 20 illustrate the spatial distributions of 
/ 7 i(r, R) at a very low energy E = 0.05 MeV, at the 
resonant energy E = 0.38 MeV, and at an energy well 
above the resonance E — 0.5 MeV, respectively. We just 
present here the distribution functions corresponding to 
the first incoming channel (K l = 0, Zj. = 0, and Z* = 0) 
which is the dominant contribution to the quadrupole 
strength function dB(E2)/dE. Other channels exhibit 
the same trends. 
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FIG. 18: (Color online) Spatial distribution of function 
/ 7 * (r, R) for E = 0.05 MeV and f = {0, 0, 0}. 
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FIG. 19: (Color online) Spatial distribution of function 
f i (r, R) for E = 0.38 MeV and f = {0, 0, 0}. 



At low energy E = 0.05 MeV, the spatial distribu- 
tion of function / 7 <(r, R) shown in Fig. 18 has a dif- 
ferent symmetry in comparison with the higher energies 
(Fig. 19, 20). We observe comparable contributions com- 
ing from two different triple-alpha configurations: the 
prolate triangle and the oblate triangle as shown in Fig. 9a 
and Fig. 9b, respectively. There is a large cancellation be- 
tween these two contributions resulting in a small value 
of dB(E2)/dE at low energies (for example, the two con- 
tributions are ~ 10" 12 but their sum is ~ 10 16 ). Never- 
theless these cancellations are well within the numerical 
accuracy of our computations. 

Fig. 19 and Fig. 20 illustrate the spatial distribution 
of function f^i(r,R) at the resonant energy E = 0.38 
MeV and higher E = 0.5 MeV. These two cases share 
the same symmetry. The dominant contribution to the 
quadrupole strength function dB(E2)/dE comes from 
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FIG. 20: (Color online) Spatial distribution of function 
Li (r, R) for E = 0.5 MeV Y = {0, 0, 0}. 



the region within the smallest contour in Fig. 19 and 
Fig. 20 that contains both maxima B and C. The max- 
imum C in the spatial distribution of function Li (r, R) 
is caused by the triple-alpha equilateral triangle config- 
uration in both the 2^ bound state and the Hoyle reso- 
nant state (see the density distribution for each state in 
Fig. 5 and Fig. 8). The contribution to the maximum B 
mostly comes from the three-alpha oblate configuration 
which only appears in the Hoyle resonant state. Even 
though the three-alpha prolate configuration dominates 
the Hoyle state's structure, it does not contribute signif- 
icantly to the strength function because the 2+ bound 
state wavefunction is zero in that region. 



V. CONCLUSION 

We have successfully produced the triple-alpha reac- 
tion rate in the low temperature regime T < 0.1 GK 
where many numerical attempts have failed before. In 
this paper, the triple-alpha is modeled as a three-body 
Borromean system. We employ the hyperspherical har- 
monics (HH) method to tackle this problem. In the low 
temperature region, the triple-alpha proceeds through 
a quadrupole transition from the + continuum to the 
2^~ bound state in 12 C. The 2^ bound state is obtained 
by solving a set of coupled channels equations in hyper- 
radius coordinates for negative energy and the condition 
that the wavefunction goes to zero at large distances. 
Naturally the same approach cannot be applied to the + 
continuum states. We combine the R-matrix expansion 
and R-matrix propagation method in the hyperspheri- 
cal harmonics basis to obtain the solutions for the con- 
tinumm states. We also implement a technique of screen- 
ing the off-diagonal Coulomb couplings to ensure numer- 
ical stability at very low energies. 

The HHR method enables us to treat the resonant and 



the non-resonant process on the same footing. We con- 
struct the triple-alpha reaction rate from the quadrupole 
transition from the + continuum states to the 2^ bound 
state of 12 C. A thorough convergence study is performed. 
Our reaction rate converges well with screening radius 
and the size of our model space. We estimate an overall 
uncertainty of a factor of 2 for our triple-alpha rate. The 
new rate agrees with NACRE above 0.07 GK. However 
we find a large enhancement of the rate at temperatures 
below 0.06 GK which marks the transition between the 
resonant and the non-resonant process. Although we ob- 
tain a significant increase of the rate at low temperatures, 
it does not drastically change the evolution of low-mass 
stars [24] as was the case for the CDCC rates [9] . Our rate 
preserves the red giant phase and is therefore consistent 
with observations [24] . The astrophysics implications are 
not presented in this paper. A separate publication with 
more details of astrophysical impacts will appear in the 
near future. 

We also explore the importance of long-range effects 
in our problem by performing four calculations in which 
only the diagonal Coulomb couplings, the full Coulomb 
couplings, and both nuclear and Coulomb interactions 
with the off-diagonal Coulomb couplings up to 30 fm and 
800 fm are included. The results emphasize the necessity 
to include Coulomb effects correctly, especially at low 
temperatures. 

The presence of very narrow two-body and three-body 
resonances in addition to the strong, long-range Coulomb 
interaction make the triple-alpha problem very challeng- 
ing. The HHR framework which is the combination of 
various methods is a new approach to overcome the well 
known difficulty of the three charged particle system. 
This method allows us to approach the very low energy 
regime where measurements are impossible without us- 
ing extrapolation, resulting in significant impacts in as- 
trophysics. In addition, it opens new opportunities in ad- 
dressing three-body low-energy reactions in others fields 
for cases where all particles have charge. 
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Appendix A: Quadrupole transition function 

In this appendix, we present a sumary of our deriva- 
tion of the three-body quadrupole transition function in 
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hyperspherical coordinates. Our problem involves only 
spin-zero particles, notation L is therefore used to refer 
to the total angular momentum and the total spin of the 
system. A transition from the initial continuum state of 
total spin L and momentum n to the final bound state 
of spin L' is characterized by the strength function: 

dB(E2, L -> V) _ E 2 /2m\ 3 
dE ~ ~2 \W) 

x f Y, \{L'M'\E 2m \LM;K)\ 2 &ni. (Al) 

The integration in Eq. (Al) is taken over all angles in 
momentum space dtt£ = sin# K cos0 K d9 K dfl kx dCl ky ■ The 
final bound state wavefunction ^ L M is expanded in HH 
coordinates as: 

^'M' =p -f XK % v v (p)^ l nO)[Yv x ®Y lly ] L , M ,, 

K ' l 'J'y 

(A2) 

while the continuum state wavefunction <3/ LM has the 
following formula: 



LM 



1 



E w;M^"(w,0^]LAf 



KIJ„ 



^^{0,)[Yi % ®Y liy ]l M . (A3) 

In Eq. (A3), k is the hyper-radial momentum and rc- 
la tcs to the three-body energy E by a relationship k = 
y/2mE /ft 2 ; K l l x l l y (Kl x l y ) represents the incoming (out- 
going) channel. We first expand the squared modulus 
sum in Eq. (Al) using Eq. (A2) and Eq. (A3). The inte- 
gration in momentum space dOg is then simplified using 
the following properties: 



/ 
/ 
/ 



dn kx Y^ mi (n k jY lLm{x (n k j = s lilL 6 mimL , (A4) 

dn k Yl mi {n ky )Y l{ym{y {n ky ) = S lyl{y 5 mym{y , (A5) 
sin 2 ^cos 2 ^d^(^(0 K ))*4^(0 K ) = S K ,, K1 . 

(A6) 



We arrive at a simpler expression for Eq. (Al) which no 
longer contains the momentum space dependence 



dB(E2, L -> L') _E 2 (2m 
dE ~ ~2 { ft 2 " 



L 2 



E 



It is straight forward to obtain the expression for the 
quadrupole operator E 2m for three alpha particles in HH 
coordinates: 



E- 



eZ 



2m 



[(psmO) 2 F 2m (x) + (pcosfl) 2 Y 2m (y)] 



(A8) 
(A9) 



Inserting Eq. (A8) into Eq. (A7) we have: 
(K'l' x l' y L'\\E 2m \\Kl x l y K%l v L) s = 
^ (K'l'JyL 1 1 1 (p sin6) 2 Y 2 (x) 1 1 KlJ y K%r y L) s 

+ ^ (K'l'J' y L' \\(pcos6) 2 Y 2 (y)\\ Kl x l y K%l y L) s . 

We denote M.E.I and M.E.2 as the first and second 
term on the r.h.s of Eq. (A9). M.E.I can be factorized 
into two terms of which one contains the hyperspherical 
variable dependence and the other is angular momentum 
dependent: 



M.E.I = 



e ^{K'l'J y \{p S m6f\Kl x l y K%l y ) 
[l' x l' y L'\\Y 2 (x)\\l x l y L). (A10) 



The first bra-ket term on the r.h.s of Eq. (A10) is cal- 
culated by taking an integral over the hyper-radial and 
hyper-angular parts of the wavefunctions: 



(K%l' y \(psmd) 2 \KlJ y K%ll) = (^Pj 



-5/4 



where we define: 



Iph,e , 
(All) 



(A12) 



h fi = J sm 2 6cos 2 6d6(/* l ! y (8)ysm 2 6 4> l K lv (8). (A13) 

The second bra-ket term on the r.h.s of Eq. (A10) can 
be explicitly calculated using angular momentum algebra 
[51] as: 



(l' x l' y L' \\Y 2 (x)\\ l x l y L) = S(l y ,l' y ) x /-LP x f x 
x(-i) 2+i «+ L + 2Z x ' L ' L 2 



L X h x l y 



5 

47r" 

b x z, i x 





.(A14) 



Performing a similar calculation for the second term on 
the r.h.s of of Eq. (A9) (M.E.2), we obtain the final result 
for a three-body quadrupole transition strength function: 



E E (K'l' x l' y L'\\E 2m \\Kl x l y K%l y L) ( 



Kl x l v K'VV 



(A7) 



where the subscript s denotes the radial-space part in the 
wavefunction expansion and L 2 = 2L + 1. 



dB(E2,L^ L') 



dE 



2ft 2 \A 



E 

KHUi 



E E ipfrfiM + hfiM) 



Kl x l v K'i 



,(A15) 
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where: 



I 



IJ.. 



h.e = I sin 2 6» cos 2 6» d6» {fe? (6))* sin 2 <f> l * ly (6), 



sin 2 6» cos 2 6> d6 {<pp{6))* cos 2 9 <j>p{6), 



Ai = 8{l y ,l' v )^-l x V x {-lf +l v 
V L 2 1 ( l' x 2 l x 



2+l'+L+2l' 



1 Y Y 

L x l x l y 







A 2 = 6{l x ,Y x ) x j-i y V y {-l) 



2 + l' x +L'+l' y +ly 



L' L 2 

Ly l y l x 



Y 2 / 

ly Z ly 
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